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ABSTRACT 

Hydrogen atmosphere white dwarfs with metal lines, so-called DAZs, require 
external accretion of material to explain the presence of weak metal line absorp- 
tion in their photospheres. The source of this material is currently unknown, 
but could come from the interstellar medium, unseen companions, or relic plan- 
etesimals from asteroid belt or Kuiper belt analogues. Accurate mid-infrared 
photometry of these white dwarfs provide additional information to solve the 
mystery of this accretion and to look for evidence of planetary systems that have 
survived post main sequence evolution. We present Spitzer IRAC photometry 
accurate to ~3% for four DAZs and one DA with circumstellar absorption lines 
in the UV. We search for excesses due to unseen companions or circumstellar dust 
disks. We use Hubble Space Telescope NICMOS imaging of these white dwarfs 
to gauge the level of background contamination to our targets as well as rule out 
common proper motion companions to WD 1620-391. All of our targets show no 
excesses due to companions >20 Mj, ruling out all but very low mass companions 
to these white dwarfs at all separations. No excesses due to circumstellar disks 
are observed, and we place limits on what types of disks may still be present. 

Subject headings: circumstellar matter-planetary systems-white dwarfs 



1. Introduction 

White dwarfs have long been used to pr obe the low mass end of the IMF to look for low 



mass stellar and brown dwarf companions (IProbst fc Oconnelll Il982l ; IZuckerman fc Becklin 
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19921 ; iFarihi et al-lbood ) . With the advent of more sensitive ground- and space-based imaging 



times Jupiter's mass around white dwarfs is now possible ( 


Ignace 


2001; 


Bur lei eh et al. 


2002 


Friedrich et al. 


2005; 


Farihi et al. 


2005; 


Debes et al. 


2005b 


a 


)■ 



Searching a subset of white dwarfs that harbor markers for substellar objects can max- 
imize the return of such a survey. Nearby hydrogen white dwarfs with metal line absorption 
(DAZs) may fit this criterion. Three hypothe ses have been put forth to e xplain the presence 



of DA Zs: interstellar matter (ISM) accreti on (IDupuis et al 



20061 ). unseen comp anion wind accretion (IZuckerman et al 



poor planetesimals (lAlcock et al.lll986l ; iDebes fc Sigurdsson 



99 



2. 



1993a b; Koester & Wilken 



2003 



2002 



and accre tion of volatile 
Jural l2003h . 



ISM accretion has a wealth of problems in predicting many aspects of DAZs such as 
the large accretion rates required for some obj ects and the distribution of these objects 
with respect to known clouds of dense material (jAannestad et al.l 11993c IZuckerman fc Reid 
19981 ; IZuckerman et al.ll2003c iKilic fc Redfieldll2007l ). The quick atmospheric settling times of 
hydrogen atmospheres imply that the white dwarfs are in close proximity with accretionary 
material. 



There are roughly 40 cool DAZs known (IZuckerman et al.l 12003c iKoester fc Wilken 



20061 ). Of them, seven have dM companions, supporting the argument that DAZs could have 



unseen companion s that place material onto the WD surface through winds (IZuckerman et al. 



20031 ; iDebea 120061 ) . In order to accrete enough material, companions must be in extremely 
close orbits (P^l day), bringing into question why these objects have yet to be discovered 
through radial velocity surveys of compact objects or observable excesses in near-IR flux. 
In most cases the reflex motio n from such objects would be easily detectable , on the or- 
der of a few to tens of km/s (IZuckerman fc Becklinl 1 1 9 92c iMaxted et al.l 120061 ) . The idea 
of the presence of unseen companions also cannot explain objects like G 29-38 and 4 other 
white dwarfs which have infrared excesses due to dust disks within their host white dwarf's 



tidal disruption radius ( 


Graham et al 


1990 


2005; 


Reach et al. 


2005b: 


Kilic et al. 


2006a| 



Patterson et al. 


1991: 


Jura 


2003; 


Becklin et al. 



The disks around G 29-38 and GD 362 show 
an amorphous silicate emission featur e at ^10/xm implying a small grain size within the 
disk and possibly warped geometries (IReach et al.ll2005bl ; I Jura et al.l 120071 ). Furthermore, 
companions > 13 M j are ruled out for a wide range of orbital separations around G 29-38 
(IDebes et al.ll2005ah . 



The invocation of cometary or asteroidal material as a method of polluting WD atmo- 
spheres was developed t o explain photosph eric absorption lines due to metals in the DAZ 
WD 0208+395 (G 74-7) (lAlcock et al.lll986l ). However, the rates predicted by these original 
studies could not satisfactorily explain the highest accretion rates inferred for some objects 
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and could not easily rep roduce the distribution of DAZs based on their effective temperatures 
( IZuckerman et al.ll2003l ). However, mixing length theory predicts a drop-off of observabil- 
ity for accret ion as a func t ion o f effective temperature which may swamp out the earlier 
prediction of lAlcock et al.l (119861 ) (lAlthaus &: Benvenutol Il998l ). Also unclear is the effect 
non-axisymmetric mass loss could hav e on the fraction of com et clouds lost by their hosts 
during post main sequence evolution (IParriott fc Alcockl Il998l ). By hypothesis, cometary 
clouds are the result of planet formation, so the long term evolution of planetary sys tems 
and their interaction with these comet clouds needs to be investigated (jTremaind Il993l ) . 



The loss of a star's outer envelope during post main sequence evolution s pecifically af- 



fects t he stability of planetary systems, and can rescue the scenario proposed by lAlcock et al. 



(119861 ) . The Hill stability criterion against close approaches for two comparable mass planets 
qualitatively describes what happens to a planetary system. The stability criterion in this 
case is A c = (ai — 0*2) /a>i = 3/i 1//3 , where a is the semi-major axis, ji is the mass ratio of the 



(Hill 


1886; 


Gladman 


1993) 



loss, companions expand their orbits in a h omologous way, increasing their orbital semi- 



major axes by a factor Mj/Mj (jJeana Il924j ). The critical separation grows as the relative 
separation of the two planets stays the same, resulting in marginally stable systems being 
tipped over the edge of stability Thi s instability can lead to orbital rearrangements, the 



ejection of one planet, and collisions ( IFord et al.l 120011 ). These three events dramatically 
change the dynamical state of the planetary system. A fraction of unstable systems will 
perturb a surviving Oort cloud or Kuiper belt analogue and send a shower of comets into 
the inner system where they tidally disrupt, cause dust disks, and slowly settle onto the WD 
surface. This modification of the comet impact model can explain the accretion rates needed 
for the highest abundances of Ca observed and the presence of infrared excesses around WDs 



(iDebes & Sigurdssonll2002l ). 



The model of IDebes fc Sigurdssonl (120021 ) can be extended to asteroidal material closer 
to the star. As the central star's mass changes, the basic resonances associated with any 
planets will change and bring fresh material into unstable orbits. The amount of pollution 
will depend on the different timescales for comets and asteroids to be perturbed toward the 
white dwarf as well as the ratio of objects in either asteroidal or cometary orbits. Asteroids 
should be perturbed relatively quickly, on timescales of 10 8 yr, while comets can take up to 
an order of magnitude longer to be perturbed. Without a more detailed model, however, it 
is hard to say which population is responsible for DAZ pollution. 

Nine DAZs have already been searched for substellar companions at intermediate orbital 
separations (10 AU < a < 50-100 AU) with NICMOS high contrast imaging and AO imaging 
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dKuchner et"aDll998l : [Debes et al.ll2005bl Jal. l2006h . No planets >10 Mj were detected for four, 
and no brown dwarfs >29 Mj were detected for the other five. Additionally, no unresolved 
companions were detected down to substell ar limits, following a general finding for a dearth 



of substellar obje cts around white dwarfs (IFarihi fe Christopher 
Farihi et al.ll2005h . 



2004; Dobbie et al. 2005 



With the launch of Spitzer an unprecedented sensitivity is now possible to further con- 
strain the presence of companions in close orbits, as well as the presence of dusty disks. A 
large interest in infrared excesses aro und white dwarfs in general is evidenced by the many 



surveys of white dwarfs with Svitzer ([Hansen et al 



20061 ; Ivon Hippel et al.l 120071 ; lJura et al 



2006 



20071 : 1 Jura et all 120071 ) 



Kilic et al.ll2006bl ; Mullallv et al. 



In this paper we present results of our search of four nearby DAZs and a DA with 
circumstellar absorption that have no known excesses for companions and circumstellar disks. 
In §2] we detail our Spitzer IRAC photometry and results, while in $3] we present second 
epoch NICMOS images of WD 1620-391 to look for common proper motion companions to 
the white dwarf. Finally in £J3]we present our conclusions. 



2. Spitzer Photometry 
2.1. Observations 

Table [T] shows our target DAZs, complete with known T e ff, log g, distances, and ages. 
Cooling ages were taken from the li terature and initial masses and main sequence lifetimes 
were calculated by the equations of IWoodl (J1992J): 
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10.4 In 



t 



MS 



M WD 
0.49M P 









lOM,(M )- 2 - 5 Gi/r. 



(1) 
(2) 



Each target was obse rved with the fou r IRAC channels, with nominal wavelengths of ~3.6, 
4.5, 5.8, and 8.0 /xm (jFazio et al.ll2004l ). The observations were carried out in the mapping 
mode, with 30 random point dithers for each pair of channels. At each dither point, the 
camera integrated for 100 s, for a total of 3000 s in each band. The exception to this was 
WD 1620-391, which is a much brighter source. The images had exposure times of 30 s per 
dither with 75 dithers for a total integration of 2250 s. Table [2] summarizes our observations. 

In order to obtain Spitzer IRAC photom etry with an accuracy of ~3%, we followed the 
prescription laid out in lReach et al.l (j2005al ). We took the BCD files from the latest Spitzer 
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pipeline calibratio ns for each target (S14.0) a nd created a final, mosaicked image using the 



MOPEX package (IMakovoz fc Marleaul 120051 ) . Some caution for point source photometry 



with IRAC is warranted. Post-BCD pipeline calibrated mosaics are not of a high enough 
fidelity for accurate photometry of stellar point sources. We routinely found that PBCD 
images returned photometry systematically 2-4% higher than when we used MOPEX. We 
performed overlap correction with a default overlap correction namelist, and mosaicking with 
the default namelist given in the IRAC data handbook. For brighter point sources, the outlier 
rejection schemes of MOPEX can spuriously reject good pixels as cosmic rays due to photon 
noise larger than the background variation. A typical symptom of this is a coverage map 
file that shows that many images were thrown out at the position of the target source. We 
experienced good results by choosing an UPPER_ and LOWER_THRESHOLD parameter of 
15 for the MOSAICIN module, as well as using the keyword REFINE_OUTLIER to ensure 
bright sources were treated with a threshold closer to 20. The thresholds refer to the number 
of sigma above the mean background. As a final check we visually inspected the resulting 
coverage maps to ensure that most images were used by the mosaicking program. 

Since each of our images had several dither positions, we did not make any array- 
location or pixel phase corrections. We estimate that these effects are at the level of 1% 
and not a significant error source, but we include them in our total error. We performed 
aperture photometry with a 3 pixel radius (~3.6"), and used a 4-pixel wide annulus starting 
just outside the source aperture for background subtraction, to ensure as accurate estimate 
of the background as possible. Aperture corrections appropriate for this size source radius 
and background annulus were applied, as well as calibration factors, flux con versions and 



a color correction in each band assuming a v spectral slope as mentioned in iReach et al. 



( j2005al ). The consistency of both aperture corrections and the photometry with different 
sized apertures was checked by recalculating the photometry with 5 pixel radius apertures 
with background annuli with 5-pixel radii starting just outside the source aperture, and 3- 
pixel source apertures with 10-pixel wide annuli starting at a radius of 10 pixels. We avoided 
a 2-pixel source aperture as that appeared to consistently give photometry lower by ~2-5%. 
For channels 1 and 2, differences between the three choices were never more than 1% except 
in the case of WD 0245+541, which has several nearby sources within 4-10 pixels. Channels 
3 and 4 often had larger changes for the 5 pixel radius aperture, up to ~10% but typically 
closer to 2%. We attribute these systematic changes primarily to residual structure in the 
background and to coincident sources. Both of these sources of systematic error are lessened 
by the small aperture and small background annulus. We estimate that on average there is 
a 1% error from sytematic uncertainties in aperture photometry based on our specific choice 
of aperture and background annulus. 

No obvious interstellar cirrus was noted for any of our targets in the 8/xm channel. 
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Figures[T]-[3]show PSF subtracted NICMOS images of the DAZs from lDebes et al.l (l2005bf ). 
with contours from the final IRAC channel 2 images overlaid. The contour lines correspond 
to 0.1%, 1%, and 10% of the total measured flux respectively to demonstrate the absence 
of contaminating objects in the source and background photometric apertures. The target 
WD in each image is located at the point (0,0), and appears as a speckled area since it is 
behind the coronagraph and the residual PSF has been subtracted off. 

For the observations of WD 0208+396, the IRAC detector was struck by a large number 
of solar protons, degrading the images with cosmic ray hits. The looser constraints on outlier 
rejection can give higher counts at the level of 10%. These hits were worse for the 5.8/xm 
channel but we used a more stringent threshold for the MOPEX outlier routines of 3 for 
channels 3 and 4 instead of 15. Inspection of the coverage maps for the channels show that 
most of the images could still be used, with the most images being rejected for the 5.8/xm 
channel. We verified that we got consistent photometry by visually inspecting individual 
BCD images and combining only the files without obvious cosmic ray strikes. 

The estimated photometric errors for each channel are quite small due to the large 
S/N achieved. In addition to the standard errors in photometry, we added a 3.3% factor to 



account for the overall uncertainty in the flux calibrations quoted by lReach et al.l (j2005al ) as 
well as the contributions from uncertainties mentioned above. 



2.2. Comparison of Photometry to WD models 



In order to detect a bona fide excess, one must com pare the observed flu x with an 
expected flux. We compared ou r observations w i th mo dels of iBergeron et al.l (119951 ) as well as 
the BVRIJHK p hotometry of IBergeron et al.l (120011 ) for four of the five targets. WD 1620- 
391 was not part of IBergeron et al.l (120011 ) 's survey and so we used a combination of USNOB, 
Hipparcos, and 2MASS photometry. Fluxes in the mid-infrare d were kindly provided (P.E . 
Tremblayprivate communication), using updated models from iTremblay fc Bergeron! (120061 ) 
and without any knowledge of the measured mid-IR fluxes. We further normalized these 
flux densities to a median of the visible and Near-IR flux densities to account for any slight 
offsets between the observed data and the models. This approach differs from previous work 
reported, where blackbody extrapolations o f the WDs' K flux density were compared with 
our Spitzer data (IDebes fc Sigurdssonll2007l ). 



For the level of photometric accuracy we have achieved, white dwarfs with effective 
temperature of ~5000-7000 K depart from true black bodies, mainly due to H _ bound- 
free and free-free opacity, with the free-free opacity being most important for the near- and 
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mid-infrared (P.E. Tremblay, private communication). Free- free abs orption can be calculated 



precisely at long wavelengths and is incorporated in WD models (see lJohnlll988l . for example 



Figure H] shows a representative comparison between the model fluxes and the measured 
fluxes for WD 0208+396, as well as the residuals. The full list of predicted and observed 
IRAC fluxes for all of our targets is in Table [3j while Figures [5] and [6] show the SEDs of 
the remaining targets. We required that a significant excess (deficit) be > three times the 
photometric error above (below) the calculated model flux in at least one channel. We find 
that for the exception of WD 1257+278, the model fluxes and photometry agree to within 
1-2 a. 

Figure [6] shows the SED of WD 1257+278 compared to the model. There isexactly a 3<r 
deficit in the 4.5 band, to a depth of 10%. The mosaic coverage maps show no images being 
thrown out where the photometric aperture is located. A slight mismatch between the model 
effective temperature and the true effective temperatrue could present an artificial deficit or 
excess, but the errors in the derived effective temperature are on the order of ~2%, which 
would correspond to errors in the predicted fluxes of 3-4%, much less than the observed 
deficit (P.E. Tremblay, private communication). Despite matching our criteria for selection 
as a significant deficit, we believe it is tentative at best, based on a detailed analysis of the 
match between our photometry and the models. 

Because of the deficit with WD 1257+278 we wished to get an empirical sense of how 
well the data matched the predicted model fluxes. To that end, we took the standard 
deviation of AF U /F U ^ P in all the channels where AF U is the difference between the observed 
flux density and the predicted flux density (F UtP ), as well as the mean AF u /F UyP for each 
channel. We find that the standard deviation of the sample is ~3.7%, while the mean for 
each channel is -l%,-5%,0.09%, and -3%. These results indicate that the predicted fluxes 
match the observed fluxes to within the absolute calibration errors we assume. We note that 
the 4.5/im channel appears to have a barely marginal (~1.4<r) mean deficit, with four of the 
five targets possessing ~ 5% or greater deficits, WD 1257+278 being one of these objects. 
WD 0208-396 is the only object with no deficit at 4.5/im. 

As another test, we divided the IRAC photometry of our target DAZs by WD 1620-391, 
the brightest WD in our sample with the highest signal-to-noise. In this case, we are limited 
by photon noise and the stability of the IRAC detectors, which is on the level of ~2%. 
We compared the relative photometry of WD 1620-391 and WD 1257+278 to the model 
fluxes in Figure [7J Within the estimated errors, the observed flux ratios match the expected 
ratios. We repeated this test with the other white dwarfs and found similar agreement. The 
consistency of the flux ratios suggests that the depression of flux at 4.5 micron may be due 
to a systematic error in the aperture correction, color correction, or calibration factors for 
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that channel. 



Observed deficits for a white dwarf may be evidence for circumstellar material raining 
down on its surface. If such a situation were confirmed at 4.5 or 8/im, we predict that 
non-LTE absorption by SiO gas may be present, with possibly some contribution from CO. 
Absorption due to fundamental and overtone rotat ional- vibrational bands of SiO and CO 
in late type stars is well known ICohen et al.l (Il992l ). The dissociation temperature of SiO 
and CO are high enough that these species could persist at the temperatures of cooler white 
dwarfs. 

The absorption could be boosted if SiO is formed above the white dwarf photosphere 
through photodissociation of Si02 (and any CO present is similarly formed through photodis- 
sociation of carbonates) from refractory dust which sublimates as it is brought down to the 
the white dwarf surface through photon drag. The resulting SiO is formed at low densities 
just above the photosphere, and is far from local thermodynamic equilibrium, with much 
larger absorption strengths than inferred from photospheric LTE. This absorption would 
show up most strongly around 4-5/im and ~8-10//m where SiO has fundamental and first 
overtone bands at 8.0 and 4.1/zm, respective ly. CO would show up primarily in the second 
channel with its fundamental band at 4.7/im (jCohen et al.lll992l ). The details of this scenario 
need to be studied further to determine the feasibility of observing absorption due to SiO or 
CO gas. 



2.3. Limits to Companions 



For IRAC, very cool substellar objects can be detected as excesses, especially due to a 
"bump" of flux for brown dwarfs and planets at ~4.5/^m. While theoretical models predict 
the 4.5/mi flux to be large, observations of c ool brown dwarfs suggest that the spectra l 



models overestimate this flux by a factor of ~2 (IGolimowski et al.ll2004t iPatten et al.ll2006l ). 



spc 
I 



In order to place upper limits on the types of unresolved companions present around our 
targets, we compared predicted IRA C fluxes for coo l brow n dwarfs and planets by convolving 
the IRAC filters with the models of iBurrows et al.l (120031 ) appropriate for the particular age 
of each target DAZ and its distance. For the 4.5/im channel we assumed that the resultant 
flux was a factor of two smaller than predicted. We then compared our 4.5/xm 3a limits to 
those models in order to determine a mass limit. These results are presented in Table [H 



In all c ases we improve the unresolved companion limits to these objects over iDebes et al. 



fl2005bf ) by a factor of 2-4. For WD 0243-026 and WD 1620-391 we rule all companions >14 
Mj objects for separations <76 and 46 AU respectively. 
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2.4. Limits to Dusty Disks 



We can determine limits to two types of dusty disks, either geometrically flat, optically 
thick disks, such as that modeled for G 29-38 or GD 362, or diffuse, optically thin disks. 
Both GD 362 and G 29-38 can be well modeled by disks not unlike Saturn's rings, w ithin 
the t i dal radius of the wh i te dwarf with an interior edge at the dus t sublimation radius (IJura 
20031 : iBecklin et allbood : I Jura et alil200~7l : Ivon Hippel et alil2007h . 



2.4-1- Optically Thick Disks 

If we ass u me an optically thick disk, the emission of the grains can be modeled following 
Adams et all Jl987h : 

s^r^ (3) 

with T function of R: 



3tt/ V r 



(4) 



This assumes that the inner radius corresponds to a dust sublimation radius of 1200 K. 
In Table HI we show the upper limits to i based on our lack of 3<r detections in our 8.0/iin 
channel data. In most cases, excess emission would have been significantly detected at 
shorter wavelengths as well. If this type of disk is present around these DAZs, the inner 
edge of the disks must be at ;> 0.4 R , or all of them are close to edge-on. We can quantify 
the probability of observing 5 systems with inclinations determined by our upper limits out 
of a random sample of disk inclinations. For any one disk, this is ~l-cosi, and for all five 
targets the probability is negligible. Most optica lly thick dust disks observed seem to have 
exterior radii of <0.6 R (Ivon Hippel et al.l 120071 ). 



Given the 10 3 -10 4 year settling timescales (tp, See Table [Q for our targets, the lack 
of a disk does no t nece ssarily imply that the DAZs cannot accrete material in this manner. 
As lHansen et al.l (120061 ) has pointed out, the timescale for removal of dust grains within the 
tidal disruption radius of a white dwarf due to Poynting- Robertson drag is short: 



T, 



PR 



Ps 



1/zm J \3g cm 



10 10 cm 



10-3L Q 



yr 



(5) 
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where s and p s are the average grain size and density respectively, and r is the distance 
from the star, ranging from ~ 10 10 — 10 11 cm. If an incoming comet or asteroid is disrupted 
and all of the material is removed before another arrives, then some fraction of the time a 
DAZ will have this type of disk and at other times it won't while still retaining a detectable 
metal line signature. The metal line will remain detectable as long as the metal settling 
time is roughly longer than the time to the next replenishing collision. Cooler dust from 
collisions may still be detectable at longer wavelengths, or slowly drift inwards from further 
away. Using Equation 0, one can determine the rough orbital separation from which dust 
would spiral in over 1 Gyr, or a typical cooling time for a white dwarf. Assuming the typical 
values in Equation [5l dust could spiral in from as far as ~20 AU. 



2.4-2. Optically Thin Disks 



If we expect an optically thin disk, we see the emission from every emitter. If one 
assumes a particular size (and therefore a particular mass) per emitter and the number of 
emitters per unit area, one can determine the total mass in an optically thin dust disk based 
on the observed flux. We focus in particular on the limit to dust between the tidal radius 
of the white dwarf and the dust sublimation radius, since this region is of most interest for 
explaining DAZ metal accretion. 

For the sake of simplicity, we assume that a constant number density of 1 /im dust 

particles reside in a flat optically thin disk between the dust sublimation radius R su b and the 

_ i 

approximate tidal disruption radius, Rudai ~ [Pwd/ Pobj) 3 R* of the DAZ, assuming a p bj=3 
g cm -3 for the parent bodies to the dust. In this case the flux is given by a modification of 
Equation [3J 



2tcs 2 cos(z) 



n(r)B u (T)rdr 



(6) 



where we have utilized the models of lLaor fc Drainei (119931 ) to calcul ate the spherical lum 



grain temperature of each dust particle given each DAZs luminosity (IBergeron et al.l 12001 



Bragaglia et al.lll995l ). For each WD we normalize n(r) such that the resultant dust disk 
spectrum returns the 3a flux limit when convolved with the IRAC 8 /im channel filter 
response. Table H] shows the resulting upper limits for dust disk mass. For WD 1620-391, its 
radius at which dust sublimates exceeds the tidal disruption radius, and so we expect no dust 
to be present in this region. Similarly hot white dwarfs would not have dusty disks around 
them like G 29-38 or GD 362. They may have gaseous disks around them, as evidenced by 



the d iscovery of a gaseous, metal-rich disk around a hot DA white dwarf (IGansicke et al. 
2006h . 
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If there are dust disks, then dust accretion could conceivably occur for longer then the 
DAZ atmospheric settling times in our sample. However, the PR drag timescale at the tidal 
disruption radius for each DAZ is < M^sk/M. This implies that accretion is not driven by 
PR drag of a present disk. 



3. NICMOS imaging 



NICMOS coronagraphic images of these five white dwarfs were presented in lDebes et al. 



(j2005bl ). with accompanying limits to companions at 1.1/im, as well as 1.6/im for WD 1620- 
391. High spatial resolution NIR images are particularly useful for discriminating against 
potential sources of background contamination which could bias the mid-IR photometry to 
spurious excesses, given the IRAC camera's spatial resolution of 1.2"/pixel. While it may 
be rare to find coincident sources that may contaminate the photometry of the target, two 
of the five targets have visual companions within 4" of the target star. 

WD 1620-391, one of the targets with a large number of visual companions, is close 
to the galactic plane. This interesting object is not technically a DAZ. It a DA with no 
optical metal absorption lines that is a large separation common proper motion companion 
to a planet bear ing star (IMavor et a 



absorption lines (IHolberg et al 



1995: 



2004) . In the U V it possesses metallic circumstellar 



Wolff et al.ll200ll ). The planet bearing star is separated 



by 5'75 (4451 AU), and is well off the field-of-view for NICMOS. Even expecting a large 
number of coin cident sources due to its galactic latitude, it possessed an overdensity over 
that expected (IDebes et al.ll2005bl ). Motivated by this overdensity, a second epoch image 
of WD 1620-391 was obtained in March 2006, two years after the first image was taken to 
search for any common proper m otion companions. T he new image was reduced following 
the basic prescription laid out in IDebes et al.l (j2005bl ). where the white dwarf was imaged 
at two separate spacecraft roll orientations and each roll image was subtracted from the 
other and combined to produce a high contrast final image. The other objects in the field 
were masked out in the image that was used as a PSF reference, since the field of view was 
moderately crowded and subtraction residuals would hamper the detection of faint sources. 

We aligned both epochs on the pixel position of WD 1620-391 and rotated the images 
so that North was in the positive vertical direction of the images, using pixel centers and 
orientations as header keywords from the STScI pipeline. We then shifted the second epoch 
image by the measured proper motion of WD 1620-391 of 97.49±3. 28 mas/yr (n cos .^=75 .52 
mas/yr) in right ascension and 0.05±1.74 mas/yr in declination (jPerryman et al.lll997l ) to 
align the shifted background stars. We measured the centroids of ~70 observed objects 
common in both fields using the IDL ASTROLIB routine GCNTRD and measured the 
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difference in centroid position from one epoch to another. With this procedure, any object 
co-moving with WD 1620-391 would have a position shift of 2.58 NICMOS pixels, or 0"19. 

Figure [8] shows the resulting differences between the measured centroids in the two 
NICMOS image epochs. The solid circle represents the 3-cr limit as empirically measured by 
the entire sample of observed sources in the field, with la being 14 mas/yr and median proper 
motions of the sample of -17 mas/yr and -8 mas/yr. There is a slight offset in the median 
change in right ascension of the group of sources from the expected zero value, though it is 
a ~lcr difference in RA. This could be because of a bulk proper motion of the background 
sources, since WD 1620-391 is at a low galactic latitude, or a sub pixel mismatch between 
the reported pixel centers of WD 1620-391. The magnitude of centroiding errors on HST 
acquisitions, however, is closer to 7 mas and is smaller than the offset seen here. In any case, 
there appears to be no co-movin g sources, thus comp letely ruling out any companions down 



to 6Mj at separations > 13 AU (IDebes et al.l 12005 



Conclusions 



We can place stringent limits on the types of disks and unresolved companions present 
for all of our targets. For two of our targets, only planetary mass objects (M<14 Mj) can be 
present at all separations, and for the rest, only very low mass brown dwarfs (M< 20 Mj) 
can be present at separations < l"or orbital separations of between 13 and 35 AU. 

The explanation that all apparently single DAZs can be caused by the winds of unseen 
companions does not fit our results unless the companions are very low mass brown dwarfs 
or high mass planets. One would expect to see large amounts of dust present if tidally 
disrupted planetesimals or ISM accretion were the source of metals for DAZs. Our targets 
show no evidence of such dust down to ~ 10 20 g if there are optically thin disks present, 
and out to separations of ~0.4 R Q if there are optically thick disks present. We effectively 
rule out optically thick disks like those seen around G 29-38 for our targets. We cannot rule 
out dust that is further away from the white dwarf and consequently much cooler. Sensitive 
studies at longer wavelengths may yet detect dust around these white dwarfs. 

Instead, optically thick dusty disks around DA Zs seem to be somewhat rare with only 5 



such known and no optically thin disks yet reporteddZuckerman fc Becklinlll9871 ; iFarihi et al. 



2005; 


Kilic et al. 


2005, 


2006a: 



thick dust can be explained for cooler DAZs by infrequent encounters with large planetesimals 
that create short lived disks that disappear quickly while still allowing detectable metal lines. 
For that reason dusty disks should primarily be around hotter DAs, whose shorter settling 
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times require a quicker replenishment of dust and thus should have long lived disks. DAs 
that are too hot vaporize material well before it is tidally disrupted. If the disks are instead 
optically thin, then weaker emission may be present, though currently undetectable. The 
upper limits for dust disk masses imply that for many DAZs the amount of material close to 
the white dwarf is sufficient to be detectable spectroscopically, but more difficult to detect 
in the mid-IR. 
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Table 1. Properties of the Target White Dwarfs 



WD 


Name 


M f 


T e// 




D 


M 4 a 


tcool+^MS 


TD 


References 






(M ) 


(K) 


(Gyr) 


(pc) 


(M ) 


(Gyr) 


log (yr) 




0208+396 


G 74-7 


0.60 


7310 


1.4 


17 


2.1 


3.2 


3.78 


1,4 


0243-026 


G 75-39 


0.70 


6820 


2.3 


21 


3.2 


2.8 


3.39 


1,4 


0245+541 


G 174-14 


0.76 


5280 


6.9 


10 


4.6 


7.2 


4.47 


1,4 


1257+278 


G 149-28 


0.58 


8540 


0.9 


34 


1.7 


3.3 


3.26 


1,4 


1620-391 


CD-38°10980 


0.66 


24406 


0.1 


12 


3.1 


0.7 




2,3 



a See Equation [T] for the calculation of M; and the WDs' total ages. 



References. — (1) 



(4) iKoester fc WilkerJljooj ) 



Bergeron et all (1200 if ) (2) iBragaglia et all (119951 ) (3) Ivan Altena et all (1200 if ) 



Table 2. Observations 



WD 


AOR Key 


Exposure Time 


Dither Points 


Date 


Start Time 






(s) 






(UT) 


0208+396 


11389184 


100 


30 


2005-01-17 


20:35:48 


0243-026 


11389440 


100 


30 


2005-01-16 


15:44:34 


0245+541 


11389696 


100 


30 


2005-02-19 


03:34:52 


1257+278 


11389952 


100 


30 


2005-06-13 


03:18:19 


1620-391 


11390208 


30 


75 


2005-03-30 


10:12:15 
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Table 3. Predicted and Observed Fluxes in /iJy 



WD 


[3.6], 


[3.6] 


[4.5], 


[4.5] 


[5.8], 


[5.8] 


[8.0], 


[8.0] o 


0208+396 


1039 


1063+35 


669 


676+22 


426 


442+16 


238 


231+11 


0243-026 


472 


479+16 


307 


294+10 


196 


198+7 


110 


102+5 


0245+541 


1333 


1305+43 


894 


848+28 


587 


583+20 


336 


332+12 


1257+278 


300 


290+10 


192 


175+6 


122 


124+5 


68 


71+4 


1620-391 


5100 


5162+170 


3204 


3050+90 


2006 


2008+67 


1097 


1050+35 



Table 4. Excess Limits 



WD 


Companion Limit 


% A 


D b 


Dust Mass c 




Mj 


% 


R 


g 


0208+396 


20 


2.9° 


0.7 


2xl0 20 


0243-026 


14 


4.7° 


0.5 


2xl0 20 


0245+541 


20 


1.6° 


0.4 


lxlO 20 


1257+278 


20 


10.7° 


0.7 


8xl0 20 


1620-391 


13 


0.1° 


5.0 





a Upper inclination limit for optically thick disk to avoid 
detection. 

b ower limit for inner radius of optically thick disk. 
c Upper mass limit of dust for optically thin disk. 
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Fig. 1.— NICMOS F110W images of WD 0208+396 (left) and WD 0243-026 (right). The 
contours are from IRAC channel two images where the levels correspond to 0.1%, 1%, and 
10% of the total observed flux from the white dwarf. 




Fig. 2.— NICMOS F110W images of WD 0245+541 (left) and WD 1257+271 (right). The 
contours are from IRAC channel two images where the levels correspond to 0.1%, 1%, and 
10% of the total observed flux from the white dwarf. 
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Fig. 3.— NICMOS F160W image of WD 1620-391. The contours are from IRAC channel 
two images where the levels correspond to 0.1%, 1%, and 10% of the total observed flux from 
the white dwarf. 



-21 - 




Fig. 4. — Comparison of observed fluxes for WD 0208+ 396 (diamonds) and predicted fluxes 
(solid line) based on the models of lBergeron et al.l (120011 ). The bottom panel shows a close-up 
of the residuals in the IRAC channels as well as the differences compared to a pure blackbody 
SED (dashed line). 
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Fig. 5.— Same asH but for WD 0243-026 (left) and WD 0245+541 (right). 




Fig. 6- 



Same asH but for WD 1257+278 (left) and WD 1620-391 (right). 
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Fig. 7.— Comparison between the measured flux ratio of WD 1257+278 to WD 1620-391 
(squares) and that predicted by white dwarf models (dashed line). WD 1257+278 shows a 
significant deficit in its absolute photometry which is not reproduced relative to WD 1620- 
391. 
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Fig. 8. — Measured proper motion in declination vs. proper motion in right ascension for 
observed sources close to WD 1620-391. The solid circle represents the 3<r scatter of the 
observed objects, while the square denotes where an object co-moving with WD 1620-391 
would lie. 



